SU(2) symmetry in a Hubbard model with spin-orbit coupling 
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We study the underlying symmetry in a spin-orbit coupled tight-binding model with Hubbard 
interaction. It is shown that, in the absence of the on-site interaction, the system possesses the 
SU(2) symmetry arising from the time reversal symmetry. The influence of the on-site interaction 
on the symmetry depends on the topology of the networks: The SU(2) symmetry is shown to be the 
spin rotation symmetry of a simply-connected lattice, so it still holds in the presence of the Hubbard 
correlation. In contrary, the on-site interaction breaks the SU(2) symmetry of a multi-connected 
lattice. 
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I. INTRODUCTION 

The spin-orbit coupling effect as an important mech- 
anism to control spin dynamics without introducing an 
external magnetic field [l|, has received much attention 
in the context of spintronics and the attempts to build 
a spin-transistor since the first proposal by Datta and 
Das in 1990 0- Two widely discussed spin-orbit cou- 
pling contributions are the Rashba and the Dresselhaus 
eff'ects [3|-l5| . Among many interesting questions the most 
important one concerns the underlying symmetry of this 
model, which reveals many far-reaching physical implica- 
tions that are not obvious at the first glance. A paradigm 
example is the SU(2) symmetry discovered by Bernevig 
et al. in a class of spin-orbit coupled models including 
the model with equal Rashba and Dresselhauss coupling 
constants and the Dresselhaus [110] model. This find- 
ing predicted that a spin precession phenomena should 
be experimentally observable @. Most of the previous 
investigations have been focused on non-interacting sys- 
tems, while less attention has been paid to the existence 
of the electron correlations arising from the Coulomb in- 
teraction. 

In this article, we study the underlying symmetry in 
a spin-orbit coupled tight-binding model, to which only 
the time reversal symmetry is required. It is shown that, 
in the absence of the interaction between electrons, the 
system possesses the SU(2) symmetry arising from the 
time reversal symmetry. Remarkably, we find that the 
influence of the on-site interaction on the symmetry de- 
pends on the topology of the networks in the following 
way. This SU(2) symmetry is shown to be the spin ro- 
tation symmetry of a simply-connected lattice, so it still 
holds for the case of nonzero on-site interaction. In con- 
trary, the on-site interaction breaks the SU(2) symmetry 
of a multi-connected lattice. Based on the exact solution 
of a ring system, our result is demonstrated explicitly. 

The article is organized as follows: in Sec. II, we intro- 
duce a general spin-orbit coupled Hubbard Hamiltonian 



with the time reversal symmetry. In Sec. Ill, we first 
construct the SU(2) operators for an on-site interaction 
free system by using the Kramers degeneracy. In Sec. 
IV we investigate the influence of the on-site correlation 
to the SU(2) symmetry. Sec. V is the conclusion and a 
short discussion. 



II. TIME REVERSAL SYMMETRY 

The Hamiltonian H is written as follows: 

Ht=J2 4^»jCj + H.c. + J2 M^clc (1) 

i=£j i 

i 

where c[ and Ci are the creation and annihilation fermion 
operators at the ith site that have two components. 
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t''^4 



(2) 



Here Ht describes the motion of free particles, while Hu 
represents the on-site interaction of opposite spin elec- 
trons. Unlike the simple Hubbard model, Ty/ \Tij\ is no 
longer the unit matrix arising from the coupling between 
momentum (and/or position) operators and spin opera- 
tors. In this work, we do not restrict the model to be in 
a certain explicit form, we only require it possessing the 
time reversal symmetry. Therefore the conclusion of this 
article is available to the Rashba and Dresselhaus types 
of spin-orbit interactions. 

The time reversal operator for a spin-i particles takes 
the form 



r = -iayR, 



(3) 



where K denotes the complex conjugation operator sat- 
isfying 
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K (c- number) = (c- number)* K, 



(4) 



and 0°^ (a — x, y, z) are the Pauli matrices. The Hubbard 
Hamiltonian H possesses the time reversal symmetry if 
/C commutes with all the matrices Ty, i.e., [T, Ty] — 
for arbitrary {i,j}. After a straightforward algebra, we 
can find that T^ should have the form 



1 ij Zij exp 



yijiH] 



■^ 



(5) 



which can be determined by the specific model. Where 
tij is a real number, 9ij is an arbitrary angle and n^ is 
a unit vector. Obviously, Tij /tij is a unitary matrix and 
represents the spin rotation operation. 

As well known, the spin operators of the whole system 
are defined as 



N 



rcJcTn 



which obey the following commutation relations 



(6) 



(7) 



where ea/37 is the Levi-Civita symbol. 

At first, we concentrate on the case of the Hamilto- 
nian of Eq. ([ij without the spin-orbit interaction. In the 
absence of the spin-orbit interaction, 9ij — 0, we have 



^,i?e„=o] 



0, 



(8) 



i.e., the system Hg..=o possesses the SU(2) spin rotation 
symmetry. The conservation of s leads to the spin in- 
version symmetry along an arbitrary direction, which is 
in accordance with the result of the time reversal symme- 
try. According to the general analysis, the spin inversion 
symmetry is broken when the spin-orbit interaction is 
switched on. However, the time reversal symmetry still 
holds for nonzero 9ij. The objective of this article aims 
at the inverse problem, namely: Can the present sys- 
tem acquire a SU(2) symmetry from the time reversal 
symmetry? We will find that it is possible under certain 
conditions. 



III. ON-SITE INTERACTION FREE CASE 



We start with the case of zero U, but nonzero dij. Due 
to the time reversal symmetry of the present system, the 
Hamiltonian can always be diagonalized in the form 



Hj 



E 

fc.A=l,2 



^kfkxfkX =^<^kfkfk 



(9) 



where A = 1,2 labels the two fold degeneracy, k labels 
the energy levels and the corresponding two component 
fermion operators 



Jk ^ yJki ' //c2 J ' •''' ^ \ f 



(10) 



The Kramers degeneracy allows us to construct the op- 
erators 



^ 1 



fc=l 



which obey the SU(2) commutation relations 



F"F' 



'^^a.^-yi 



(11) 
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In this sense, Eq. ^ describes a non-interacting Fermi 
gas of spin-1/2 particles. The two Kramers degenerate 
states can service as the two components of a pseudo 
spin. Obviously, we have 



f,Hj 



= 0, 



(13) 



which means that the system Ht possesses a new type of 
SU(2) symmetry. Note that the construction of {F"} is 

not unique, since any linear transformation of I /^j^, f^^ ) 

cannot change the the facts of Eqs. (fT2| and (fT3|l . In 
other words, Ht is invariant under the local rotation of 
the pseudo spins {Ffe}. 

An interesting question is that: Among all the sets of 
the conserved quantities {F"}, which one closes to the 
familiar physical quantity and is feasible to measure in 
experiment. It will be shown that the geometry and the 
spin-dependent interaction play an important role in this 
issue. It is also the main goal of the present article. To 
this end, we firstly investigate the Hamiltonian Ht on 
a simply-connected network. It can be observed that, 
taking an arbitrary node as a start point there exists a 
unique path to another. A schematic illustration of the 
simply-connected network is presented in Fig. [ija). This 
characteristic feature of such networks allows Ht to be 
rewritten as 



HT^^UJdldj+Il.c. + Y,^^^dUi, (14) 

ijtj i 

by absorbing the unitary matrices in Tij into the fermion 

operators d] and dj . For instance, one can take the trans- 
formation 

^i^ij^j ~ '^ij^i^ '^ '^ ^j ~ ''ij'^i'^jj (,1^7 



by the definition 



d = c , d, = e' 



■^, 



(16) 



The equivalent Hamiltonian (|14p represents a non-spin- 
orbit interaction system. Accordingly, one can construct 
the corresponding SU(2) operators 






(17) 



satisfying the following commutation 
Similarly with Eq. (jH)), we have 



^, 



Hj 



(18) 



(19) 



The physics of o can be understood by the following 
relationship between operators iS" and sf 



S? = 



UiS^ Uj , 



(20) 



where Ui is a unitary matrix of the form e* ' ' appeared 
in the transformation 

dj = Uid. (21) 

This indicates that the operators S" act as the real spin 
operators under the local transformation {ui}. In this 
sense, one can apply the theory of itinerant electron mag- 
netism on the spin-orbit coupling system on a simply- 
connected lattice. The ferromagnetic state (all the spins 
being aligned parallel) with respect to the operators S" 
in a system with nonzero 9ij is equivalent to the spin 
helix state with respect to the operators s" Q. We will 
give an extensive discussion about this issue after taking 
the on-site interaction into account in the next section. 
On the other hand, the equivalent Hamiltonian p^ 
can be diagonalized in the form 



Ht — 2_^ ^kdi^^dka, 



(22) 



k.i7 



where 



dka — 7 ,'^jdja, 

3 

3 



(23) 
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FIG. 1: (Color online) Schematic illustration of (a) simply- 
connected network and (b) multi-connected network. If the 
connection between i and j is broken, the topology is changed 
from multi- to simply-connected one, which will affects the 
SU(2) symmetry of the system, especially in the presence of 
the on-site correlation. 



system is essentially spin rotational symmetry. The con- 
servative quantity h is connected to an experimental 
observable o , which means a persistent spin helix. It 
is hardly observed in practice since a natural material 
with simply-connected geometry is rare. However, artifi- 
cial lattices, such as arrays of quantum dots in semicon- 
ductor heterostructures [7|, |8| confining the conduction 
electrons, or optical lattices — stable periodic arrays of 
potentials created by standing waves of laser light [y], 
can implement this task. 

In contrary, for a multi-connected system illustrated 
in Fig. [ijb), the above analysis is invalid since one can- 
not find a set of unitary matrices {ui} to implement the 
transformation of Eq. (I^IT) . Then for an on-site free sys- 
tem with the time reversal symmetry, it always possesses 
a SU(2) symmetry, but the physics of the symmetry de- 
pends on the underlying topology of the network. 



Comparing Eqs. ^ and (P^ . we find that one can have 
dk = fk. (24) 

From Eq. (P^ . we obtain the identity 



^dl^dka' ^^d]^dj^,, 



which leads to 



^^f. 



(25) 



(26) 



Then we can conclude that the SU(2) symmetry obtained 
from the time reversal symmetry of a simply-connected 



IV. ON-SITE INTERACTION EFFECT ON THE 

SYMMETRY 

In previous we have found that there is a SU(2) sym- 
metry in a spin-orbit coupling system obeys the time re- 
versal symmetry when the electron-electron interaction is 
absent. This indicates the macroscopic emergence of cer- 
tain physical features (as ferromagnet, antiferromagnet, 
spin helix, etc.) in long time scale, especially in a multi- 
particle system. Nevertheless, the spin-dependent inter- 
action between particles may break the symmetry. Now 
we turn to investigate the influence of the on-site interac- 
tion Hjj to the SU(2) symmetry. Actually, applying the 
transformation di = UiCi to Hjj on a simply-connected 



system, we have 

Hu^Y. ^^c^V^tcliC^ = Yl U^4^d,^4^d,i^, (27) 



which means it is invariant under the transformation. 
Consequently for a simply-connected system, we still 
have 




(«) 



^, 



Hr 



^. 



H 



= 0. 



(28) 



This has many implications on a spin-orbit coupling sys- 
tem. Mathematically, it can be treated as a normal Hub- 
bard model. Then all the conclusions for the Hubbard 
model on a simply-connected lattice are completely avail- 
able for the present system. Here we only give a subtle 
conclusion from the Lieb theorem [10[. In the following, 
we present a statement for a Hubbard model with the 
spin-orbit interaction on a simply-connected network by 
simply modifying the abstract in Ref. [10]. In the at- 
tractive Hubbard Model (and some extended versions of 
it), the ground state is proved to have spin angular mo- 
mentum 5 = for every (even) electron filling. In the 
repulsive case, with a bipartite lattice and a half-filled 
band, the ground state has S = l/2(|i?| — |A|), where 
\B\ {\A\) is the number of sites in the B (A) sublattice. 
In both cases the ground state is unique. We believe 
that such kind of rigorous result obtained from the simple 
Hubbard model is useful for understanding the feature of 
the present model, thereby providing a general guiding 
principle for spintronics. 

Now we consider the case of a multi-connected sys- 
tem. Since the transformation di = UiCi cannot elimi- 
nate the nonzero Oij term completely, we cannot judge 



the commutation relation 



^,^1 = \t 



, Ht, 



m a gen- 



eral manner. However, a single example can provide the 
conclusion that the on-site interaction breaks the SU(2) 
symmetry generated by the operators {F"}, although the 
time reversal symmetry still holds. 

We exemplify the above analysis by taking a simple 
multi-connected network, a ring system as an example. 
This may shed light on the role of the topology of the 
network. The Hamiltonian of the ring reads 



H' 



Hi 



m\ 



rrimg . Turing 

Tlrp -t- TljJ , 



(29) 



= -J ^ c}^Cj+i - Jc\e * 2" Cat -I- H.c, 



N 



Here the spin-orbit interaction only exerts on the tun- 
neling between sites 1 and N . Despite its simplicity, it 
reveals the common properties of the underlying symme- 
try for more complex systems. It acts as a Mobius system 
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FIG. 2: (Color online) Schematic illustration of (a) spin-^ 
fermionic system of an A^'-site ring with spin-orbit interaction 
on the bond between 1 and A'', and (b) its equivalent sys- 
tem, which is a spinless fermionic system of a 2Ai'-site ring 
penetrated by a half magnetic flux quanta. 



llj . To obtain the solution, taking the transformation 



e 2« cjt 

(JV+j-l)7r 



o,. 



(30) 



e ■ 2« Cjm, = aN+j, j e [1,A^] 
the original Hamiltonian is mapped into 



rin 



2N 



e^-^^ayaj+i -f H.c, 



(31) 



It 



where we take the boundary condition o^n+i ~ ' 
represents a 2A'^-site ring penetrated by a half magnetic 
flux quanta, which is schematically illustrated in Fig. [2] 
Such a Hamiltonian can be diagonalized in the form 



njring 

rirp 



-2J^cos(/c+^)4a,, (32) 



by using the following Fourier transformation 



(33) 



where k = Tm/N, n G [1,2A'^] denotes the momentum. 
Note that the momentum shift by tt/ (2iV) in the dis- 
persion relation ensures the Kramers degeneracy. Then 
the corresponding SU(2) operators ([TT]) and the on-site 
interaction H™^ can be expressed explicitly as 



k 

k = mi/N, ne [0,N- 1], 



a-k 

0-2TT-k--n/N 



,(34) 



and 



N 






(35) 



respectively. After a lengthy but straightforward algebra, 
we have 



[r",ii""s] ^0. 



(36) 



It is clear that the time reversal symmetry is not the suf- 
ficient condition for a SU(2) symmetry, since the validity 
of the symmetry depends on the geometry of the system 
as well as the on-site correlation. Only the coexistence of 
the closed loop in the network and the on-site interaction 
between particles can break the SU(2) symmetry. 



V. CONCLUSION 

In conclusion, we studied the underlying symmetry for 
a spin-orbit coupled tight-binding model with the time 
reversal symmetry. We found that the characteristics of 



the symmetry strongly depend on the topology of the net- 
work and the on-site interaction. It is shown that, in the 
case of zero on-site interaction, the system possesses the 
SU(2) symmetry arising from the Kramers degeneracy. 
The influence of the on-site interaction on the symmetry 
depends on the geometry of the networks: The SU(2) 
symmetry is shown to be the spin rotation symmetry of 
a simply-connected lattice, so it still holds for the case 
of nonzero U. We also investigate the multi-connected 
system based on the exact solution of a simple ring. Our 
result showed that the on-site interaction can break the 
SU(2) symmetry of a multi-connected lattice. 

Regarding the reason why the topology and on-site 
correlation affects the symmetry, it may do to its gauge 
characteristic. It has been pointed that one can regard 
the Rashba and the Dresselhaus spin-orbit interaction 
in two-dimensional semiconductor heterojunctions as a 
non-Abelian gauge field, or the Yang-Mills field U^i . The 
Yang-Mills field generates a physical field due to which 
the wave function acquires a spin-dependent phase fac- 
tor. Therefore it is not surprising that the topology of 
the system affects the symmetry. 
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